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INTRODUCTION 


Composite  solid  propellants  consist  essentially  of  oxidizer  particles  of 
varying  sizes  interspersed  in  a  binder  matrix.  The  combustion  of  such  a  system  is 
phenomenologically  complex.  A  number  of  physico-chemical  processes  are 
known  to  occur.  However,  detailed  knowledge  of  these  processes  is  not  currently 
available.  Very  significantly,  many  of  these  processes  occur  in  a  spatio-temporal 
extent  which  virtually  forbids  any  meaningful  experimental  observation  at  the 
combustion  mechanisms.  In  order  to  achieve  any  degree  of  predictive  capability, 
analytical  modeling  of  this  system  is  necessary.  A  rigorous  mathematical 
description  of  the  combustion  of  composite  propellants  represents  a  formidable 
problem. 

In  an  effort  to  circumvent  some  of  the  difficulties  inherent  in  observation 
or  analytical  modeling  of  the  combustion  zone  produced  by  the  chaotic 
microstructure  of  a  typical  composite  propellant,  studies  have  been  made  at  edge 
burning  laminates.  This  configuration,  called  a  "sandwich,"  consists  of  a  layer  of 
hydrocarbon  fuel  between  two  layers  of  oxidizer.  This  permits  at  least  some 
fundamental  aspects  of  composite  propellant  combustion  to  be  reproduced  within 
the  framework  of  a  simpler  geometry  that  is  easier  to  model  analytically,  and 
easier  to  observe  experimentally  than  is  possible  with  a  propellant. 

There  is  a  large  body  of  work  on  sandwich  combustion.  There  are  several 
experimental  investigations  (Ret.  1-14).  These  investigations  include 
cinephotographic  observations  of  burning  sandwiches,  optical/scanning  electron 
microscope  observations  of  quenched  samples  and  schlieren  observation  of  the 
gas  phase  structure.  Several  results  are  reported  from  burning  rate 
measurements  tor  different  oxidizer-fuel  combinations  and  laminate  thicknesses 
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(Ref.  2,3,5-8,11).  The  combustion  of  sandwiches  has  been  studied  at  several 
pressures. 

While  it  is  not  clear  that  the  microscopic  flame  is  revealed  in  optical 
observations  of  combustion,  its  nature  is  suggested  by  visible  density  fluctuations 
above  the  flame  (Ref.  6,10-12).  High  speed  photography  shows  that  this  is  related 
to  a  flickering  or  local  intermittency  along  the  flame  sheet  (Ref.  11).  At 
pressures  below  the  AP  self -deflagration  limit  (  —  2  MPa),  the  flow  above  the 
flame  is  laminar.  Systematic  investigation  of  the  surface  profile  of  burning 
sandwiches  has  been  made  via  the  interrupted  burn  technique.  The  leading  edge 
of  regression  is  found  to  be  in  the  oxidizer  portion  of  the  sample  under  several 
conditions  (Ref.  4,8,10,13).  A  singular  surface  condition  characterized  as  the 
"smooth  band"  has  been  observed  in  the  oxidizer  region  close  to  the  oxidizer- 
binder  interfaces  (Ref.  13,14).  Sandwich  deflagration  limits  have  been  found  (Ref. 
13).  This  is  characterized  as  the  binder  thickness  below  which  the  sandwich  does 
not  bum  at  a  given  pressure. 

There  are  several  analytical  investigations  of  sandwich  combustion  (Ref. 
15-19).  A  comparison  of  the  analytical  and  experimental  investigations  of 
sandwich  combustion  reveal  certain  gaps  of  analytical  modeling  to  date.  The 
analytical  approaches  have  not  succeeded  in  integrating  into  them  even  the 
widely  recognized  features  of  sandwich  combustion.  In  particular,  adequate  use 
of  the  information  on  the  surface  profile  geometry  has  not  been  made.  A  clear 
appreciation  does  not  exist  of  the  length  scales  involved  in  the  dynamics  of 
sandwich  combustion.  Most  representations  of  the  flame  structure  used  in 
modeling  are  simple  geometric  entities  which  are  unrealistic  in  view  of  the 
experimental  evidence.  Consequently,  even  the  basic  features  of  the 
microstructure  of  the  flame  complex  remain  unrecognized  by  theoreticians  and 


unmeasured  by  experimentalists.  No  quantitative  theory  exists  that  relates  the 
structure  of  the  flame  complex  and  the  details  of  the  surface  profile.  Generally, 
simplistic  representations  are  made  of  the  fluid  dynamical  problem  which  is  an 
inevitable  concomitant  of  the  thermal  problem.  It  follows  that  no  realistic 
estimates  exist  of  the  magnitude  of  the  convective  heat  transfer  in  the  problem. 
The  present  work  attempts  to  address  these  issues. 

It  is  easy  to  recognize  that  a  model  which  incorporates  all  the  features  of 
sandwich  combustion  would  be  a  very  complex  entity.  The  solution  of  such  a 
system  is  held  to  be  very  difficult.  The  current  model  considers  the  simpler  two- 
slab  geometry.  The  philosophy  of  the  current  approach  was  to  include  as  many  of 
the  experimentally  observed  features  as  deemed  consistent  with  mathematical 
tractability.  It  was  considered  particularly  important  to  include  information  on 
the  surface  profile,  a  quantity  which  can  be  observed  with  greater  spatial 
resolution  than  any  other  relevant  feature  of  the  combustion  zone  (because  of 
quench  testing).  The  present  approach  is  an  inverse  to  the  direct  modeling  of  the 
combustion,  in  that  the  experimentally  determined  surface  profile  and  sample 
burning  rate  are  used  as  inputs  to  a  calculation  aimed  at  determining  features  of 
the  flame  complex.  The  analytical  model  was  primarily  motivated  by  the  need  to 
recover  quantitative  information  about  the  gas  flame  microstructure.  The 
specific  objects  of  interest  in  the  investigation  are: 

1)  the  characteristic  length  scales  of  the  problem, 

2)  the  nature  of  the  flame  structure,  i.e.,  the  magnitude  and  distribution  of 

the  heat  release  rate, 

3)  the  role  of  fluid  dynamics  of  the  configuration,  and 

4)  the  sensitivity  of  the  flame  structure  to  a)  changes  in  the  surface 

profile,  and  b)  the  physical/  chemical  properties  that  characterize  the 


system. 


Experimental  work  was  also  initiated,  its  object  being  the  measurement  of 
the  surface  profile  of  the  oxidizer  and  the  fuel  and  the  interface  burn  rate. 

MODEL  DESCRIPTION 

The  model  is  concerned  with  the  combustion  at  the  interface  of  a  semi¬ 
infinite  slab  of  oxidizer  and  one  of  fuel.  The  shape  of  the  typical  burning  surface 
profile  is  shown  in  Fig.  1.  This  figure  is  representative  of  the  geometry 
considered  in  the  model.  The  origin  of  the  coordinate  system  is  at  the  tip  of  the 
cavity  and  it  rides  the  surface  at  the  interface  burning  rate.  The  functions  ysj(x) 
and  ys2(x)  describe  the  surfaces  of  the  oxidizer  and  the  fuel  respectively,  and  n  is 
the  unit  normal  vector  (positive)  pointing  outward  into  the  gas  phase. 

The  principal  assumptions  of  the  model  are  as  follows: 

1)  The  system  achieves  steady  state.  Under  this  condition,  the  geometry  of 
the  surface  profile  is  not  time  dependent.  There  is  experimental  evidence  to 
substantiate  this  claim  (Ret .  6,10-12). 

2)  The  model  assumes  that  surface  melts  have  negligible  effect.  This 
assumption  restricts  the  interpretation  of  the  results  from  the  model  to  oxidizer- 
fuel  systems  that  exhibit  no  or  minimal  melts  under  the  conditions  modeled. 

3)  The  oxidizer  and  the  fuel  are  assumed  to  pyrolyze  at  the  surface  in  a 
manner  describable  by  the  classical  Arrhenius  rate  law. 

4)  The  oxidizer  is  assumed  to  be  a  non-monopropellant.  The  model  would, 
however,  be  applicable  to  AP  burning  in  this  system  below  its  deflagration  limit. 
This  restriction  is  only  to  ensure  that  the  system  burns  in  a  geometry  typified  by 
the  one  in  Fig.  1,  used  in  the  computational  program. 

5)  All  processes  occurring  in  the  system  are  two-dimensional.  Several 
additional  assumptions  have  been  made  to  ensure  computational  simplicity. 


These  are  discussed  in  the  context  of  their  application.  The  equations  and  the 
boundary  conditions  which  describe  the  system  are  discussed  next. 

Equation  (l)  is  the  two-dimensional  energy  equation  for  the  gas  phase. 
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It  is  assumed  that  Le  =  1  everywhere  in  the  gas  phase.  Heat  transfer  is  assumed 
to  occur  due  to  temperature  gradients  only.  The  variation  of  gas  density  with 
temperature  is  considered.  The  specific  heat  capacity  and  the  thermal 
conductivity  of  the  gas  mixture  (as  well  as  those  of  the  solid  oxidizer  and  the 
solid  fuel  in  Eq.  (2)  and  (3))  are  considered  to  be  constants.  The  right  hand  side 
forcing  function  in  Eq.  (1)  is  the  heat  release  rate  per  unit  volume.  It  is  the 
unknown  which  is  sought.  It  may  be  pointed  out  that  the  description  of  the  gas 
phase  chemical  reaction  rate  is  implicitly  lumped  into  this  term.  Individual 
species  transport  equations  are  not  solved.  It  was  thought  that  primary  emphasis 
of  this  investigation  should  be  directed  toward  energy  distribution,  rather  than 
the  chemistry  of  the  process. 

Equations  (2)  and  (3)  are  the  energy  transport  equations  for  the  solid 
oxidizer  and  the  fuel,  respectively. 
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It  is  a  requirement  tor  steady  state  that  there  be  no  lateral  velocities  in  the 
solids,  and  hence  no  term  for  convection  in  the  x  direction  in  Eq.  (3). 

Equation  (4)  is  the  statement  at  mass  conservation  tor  the  two-dimensional 

tlow. 

3(pu)  +  9(pv>  =  o  (4) 

Ox  Oy 

The  variation  at  density  is  accounted  tor.  The  density  varies  not  due  to 
compressibility,  but  due  to  variation  ot  temperature  in  the  field.  The  tlow  is 
incompressible. 

Equations  (5)  are  the  Navier-Stokes  equations  tor  the  tlow.  It  has  been 
assumed  that  the  tlow  is  laminar,  an  assumption  that  is  supported  by  the  low 
Reynolds  number  and  by  experimental  results  tor  sandwiches  burning  at  low 
pressures  (Ret.  6,10-12). 
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The  velocity  and  the  temperature  fields  are  coupled  by  virtue  of  the  variation  of 
density  with  temperature.  This  has  very  important  consequences  for  the  entire 


model. 


Equations  (6)  are  the  classical  Arrhenius  rate  laws,  which  describe  the 
pyrolysis  taking  place  at  the  surf  ace  of  the  oxidizer  and  the  f  uel. 
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The  constants  in  these  equations  are  treated  as  known  inputs  to  the  model. 


Equations  (7)  express  the  appropriate  tar  tield  boundary  conditions  in  the 


condensed  phase. 
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Equation  (8)  is  a  representation  of  the  constancy  at  the  heat  flux  on  either 
side  of  the  solid  oxidizer  -  solid  fuel  interface. 
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tor  x  >  0  ;  at  the  surface 


It  has  been  assumed  that  radiant  transport  at  energy  is  not  important.  It  has  also 
been  assumed  that  the  kinetic  energy  at  the  interface  is  negligible.  The  surface 
processes  that  occur  may  either  be  endothermic  or  exothermic. 

Equation  (10)  is  a  statement  ot  mass  conservation  at  the  above  interfaces. 


/  A  A 


A  .  A 


pg<u»  = 


Pj  r  ny  tor  x  <0 


P2  r  ny  tor  x  >0 


Equations  (11)  describe  the  geometry  ot  the  problem. 
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It  may  be  noted  that  a  non-symmetrical  protile  has  been  provided  tor. 


APPROACH  TO  THE  SOLUTION 


The  system  of  equations  dealt  with  are  elliptic  partial  differential 
equations.  Equation  (1)  is  nonlinear  by  virtue  of  the  variation  of  gas  phase 
density  with  temperature.  Equations  (5)  are  nonlinear.  The  domain  on  which  the 
equations  are  valid  is  of  rather  complicated  geometry.  Based  on  these 
observations,  it  was  decided  to  adopt  a  numerical  approach. 

At  the  outset,  the  surface  profile  geometry  h  (x,y)  and  the  interface  burn 
rate  are  specified.  The  pressure  at  which  a  steady  state  configuration 
characterized  by  the  above  is  also  known.  The  observation  is  made  that  the  gas 
velocity  vectors  at  the  solid-gas  interface  are  locally  normal  to  the  wall.  This  is 
a  requirement  of  no  slip  flow.  Using  the  above  and  Eq.  (10),  the  velocity  vector 
components  at  the  wall  are  determined.  If  the  temperature  field  is  assumed,  the 
density  field  can  be  determined.  Upon  a  suitably  defined  closed  domain  (with 
appropriately  defined  bondary  conditions)  the  Navier-Stokes  equations  can  be 
solved  to  obtain  the  velocity  field  in  the  gas  phase. 

The  specification  of  the  interface  burn  rate  and  the  surface  profile  results 

A 

A  A 

in  the  knowledge  of  the  surface  temperature  distribution  T$  (x,y),  as  seen  from 
Eq.  (6).  This,  in  turn,  permits  the  solution  of  the  energy  equations  (2)  and  (3)  in 
the  solid  phase.  The  boundary  conditions  (7)  and  (8)  are  used  in  this  process.  The 
condensed  phase  temperature  field  is  thus  computed.  The  computation  of  the 
temperature  gradients  at  the  surface  in  the  solids  is  now  possible. 

We  turn  next  to  the  gas  phase  energy  equation,  Eq.  (1).  The  velocity  field 
is  known  as  a  result  of  the  solution  of  the  Navier-Stokes  equations.  A  plausible 
heat  release  rate  distribution  is  assumed.  The  surface  temperature  distribution 
has  already  been  determined.  This,  along  with  other  boundary  conditions,  leads 


to  a  solution  of  the  temperature  field  in  the  gas  phase.  The  density  field  can  now 
be  computed  as  well  as  the  temperature  gradients  at  the  surface  in  the  gas.  Use 
of  Eq.  (9)  is  now  invoked.  The  satisfaction  of  this  condition  to  a  specified 
accuracy  is  used  as  a  matching  condition.  The  process  of  solving  Eq.  (1)  is 
repeated  using  improved  estimates  of  F(x,y)  until  the  required  match  is  obtained. 

It  might  be  pointed  out  that  a  further  overall  iterative  process  involving 
the  solutions  to  the  gas  momentum  and  energy  equations  is  necessary.  This  is 
due  to  coupling  between  the  two  systems  via  the  gas  phase  density.  An 
appropriate  overall  convergence  criterion  may  be  defined  f or  this  process. 

NUMERICAL  TECHNIQUES 


Navier-Stokes  Equations 

Experimental  results  point  to  the  solid  oxidizer  and  the  solid  f  uel  receding 
less  and  less  rapidly  farther  away  from  the  tip  of  the  cavity.  It  has  been  assumed 
that  the  profile  approximates  a  channel  geometry  far  from  the  tip  of  the  cavity. 
The  geometry  considered  for  the  fluid  problem  is  shown  in  Fig.  2.  The  walls  BC 
and  DE  are  no  slip  boundaries.  BAE  is  the  combustion  driven  inflow  boundary.  A 
Poiseuille  boundary  condition  is  initially  used  tor  the  outflow  boundary  (though 
subsequent  updates  tor  this  boundary  are  based  on  a  "tree"  condition).  The 
location  at  which  this  boundary  condition  is  applied  is  decided  by  numerical 
experiments. 

The  boundary  is  transformed  into  a  rectangle  using  the  non-coni ormal 


transformation: 


ysl  (x)  tor  x  >  0 
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y$2  (x)  tor  x  <  0 


The  equations  are  non-dimensionalized  in  the  conventional  manner,  and  the 
problem  is  recast  in  vorticity-stream  function  variables.  The  solution  can  then 
be  obtained  tor  a  known  density  distribution  in  the  field.  The  nonlinear  vorticity 
transport  equation  is  solved  by  successive  line  under  relaxation  technique  with 
linear  extrapolation  tor  the  nonlinear  terms.  The  boundary  vorticities  are 
updated  by  first  order  representations.  The  Poisson  equation  is  solved  by  the 

A 

successive  line  over  relaxation  method.  A  variable  grid  is  provided  in  the  n 
direction.  The  solutions  to  the  vorticity  transport  equation  and  the  Poisson 
equation  are  successively  iterated  within  an  overall  iteration.  The  details  of  all 
computations  in  this  work  are  reported  elsewhere  (Ret.  20).  The  solutions  were 
obtained  on  a  CDC  730/760  system.  Typical  computation  times  were  on  the 
order  of  150  CPU  secs.  The  streamlines  were  plotted  by  an  auxiliary  code. 

Condensed  Phase  Energy  Equations 

The  solution  to  the  energy  equation  in  the  solid  oxidizer  and  the  solid  fuel 

A  * 

are  obtained  in  tandem.  The  length  scales  are  non-dimensionalized  by  a^/r  and 

A 

the  temperature  by  Tq.  The  equations  are  then  transformed  using  Eq.  (12).  The 
geometry  and  the  boundary  conditions  for  these  equations  are  shown  in  Fig.  3. 
The  numerical  boundary  conditions  corresponding  to  Eqs.  (7)  are  applied  at  finite 
distances  whose  locations  are  decided  by  numerical  experiments.  A  special 
treatment  of  the  code  is  necessary  for  the  vertical  interface.  The  solution  is 
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Figure  3.  Geometry  and  boundary  conditions  for  the  condensed 
phase  energy  equation. 
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obtained  by  successive  line  over  relaxation.  The  temperature  tield  obtained  is 
displayed  as  a  plot  of  isotherms.  Typical  computational  time  on  the  CDC 
760/730  system  is  about  40  CPU  secs. 

Gas  Phase  Energy  Equation 

The  density  in  Eq.  (1)  is  replaced  by  pressure  and  temperature  by  using  the 
perfect  gas  law.  The  equation  is  then  non-dimensionalized  and  transformed  as 
before.  The  transformed  domain  is  a  rectangle  as  for  the  fluid  problem.  The 
equation  can  then  be  solved  if  the  velocity  tield  is  known  (from  the  solution  of 
the  Navier  Stokes  equations)  and  a  heat  release  distribution  is  assumed.  The 
nonlinear  terms  are  linearized  by  extrapolating  the  values  from  previous 
iteration  levels.  The  equation  is  solved  by  the  successive  line  under  relaxation 
method.  The  temperature  field  is  displayed  as  a  plot  of  isotherms.  A  typical 
solution  takes  about  100  CPU  secs  on  the  CDC  760/730  system. 

MODULAR  RESULTS 

The  baseline  objective  of  this  work  is  to  obtain  information  on  the  gas 
phase  combustion  zone  microstructure.  This  involves  the  solution  of  several 
equations  in  conjunction  with  each  other  within  an  overall  iterative  process.  This 
procedure  is  not  conducive  to  obtaining  insight  into  the  behavior  of  particular 
aspects  ot  the  system.  It  was  therefore  decided  to  work  first  with  smaller,  more 
easily  comprehended  parts.  The  condensed  phase  energy  problem,  the  gas  phase 
momentum  problem  and  the  gas  phase  energy  problem  constitute  such  modular 
entities.  The  object  is,  in  each  case,  to  delineate  a  general  set  ot  rules  ot 
behavior  valid  within  a  framework  of  constraints.  It  is  acknowledged  that  care 
must  be  taken  in  interpreting  these  results  in  order  that  no  more  significance  be 
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attributed  to  them  than  is  warranted  under  the  imposed  constraints.  Since  the 


analytical  work  was  initiated  prior  to  the  experimental  part  of  this  investigation, 
the  surface  profile  and  the  interface  burn  rate  from  the  investigation  of 
Bakhman  and  Librovich  (Ret.  7)  corresponding  to  the  combustion  of  KP-PMMA 
slab  system  at  a  pressure  of  30  atm.  was  used  tor  the  f  ollowing  modular  results. 


Fluid  Module 

The  code  was  initially  run  with  constant  density  in  the  entire  field.  The 
velocity  field  so  obtained  was  used  in  the  computation  at  the  solution  of  the  gas 
phase  energy  equation  along  with  an  assumed  nominal  heat  release  distribution. 


From  this  the  variable  density  field  was  computed.  The  solution  to  the 
momentum  equations  was  again  recomputed  for  the  variable  density  case.  The 
Re  tor  this  case  was  58.  Figure  4  shows  the  results  obtained  for  the  case  with 
constant  density  as  well  as  the  one  with  variable  density.  It  is  a  plot  of 
streamlines.  For  the  variable  density  case  the  flow  preserves  its  two- 
dimensional  character  farther  out.  Also,  the  curvature  of  streamlines  passing 
through  the  flame  region  is  discernible.  Apart  from  these  qualitative  changes, 
there  are  significant  changes  in  the  magnitudes  of  the  velocities  for  the  two 
cases.  The  coupling  between  the  momentum  and  the  energy  equations  is  thus 
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strong. 

Condensed  Phase  Module 

For  representative  values  of  thermophysical  properties  of  KP  and  PMMA 
the  solution  obtained  for  the  condensed  phase  energy  equation  is  shown  in  Fig.  5. 
The  temperatures  in  the  figure  are  nondimensionalized.  In  the  oxidizer  the 
temperature  gradients  in  the  x-direction  are  considerably  less  than  those  in  the 
y-direction.  The  isotherms,  in  general,  tend  to  follow  the  surface  profile.  The 


Streamline  plot  of  the  baseline  fluid  module  solution, 


Baseline  condensed  phase  isotherm  plot. 
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direction  ot  heat  transter  at  the  vertical  interlace  is  toward  the  luel  at  distances 
at  the  order  of  50  |im  from  the  tip  of  the  cavity.  At  a  distance  of  about  200  /in1 
the  vector,  though  very  small  in  magnitude,  is  pointed  toward  the  oxidizer.  The 
effect  ot  thermophysical  properties  of  the  oxidizer  and  the  luel  has  been  studied 
in  detail  (Ref.  20).  These  aflect  the  solution  significantly  in  respect  to  the  above 
considerations.  The  effect  of  the  surface  profile  was  also  studied.  It  has  been 
concluded  that  the  slopes  of  the  fuel  and  the  oxidizer  at  the  tip  of  the  cavity 
significantly  affect  the  magnitude  and  the  direction  of  the  heat  flux  vector  at 
the  vertical  interface.  The  effect  of  pyrolysis  kinetics  on  the  temperature  field 
enters  via  the  surface  temperature  distribution.  In  general,  the  surface  profile 
effects  are  more  significant  than  the  surface  temperature  distribution  effects. 

Gets  Phase  Module 

Figure  6  shows  the  baseline  result  obtained  tor  the  gas  phase  energy 
equation  with  a  heat  release  distribution  assumed  tor  a  nominal  match  of  the 
temperature  gradients  at  the  surface.  It  is  a  plot  at  isotherms.  The  results  point 
to  the  bulk  of  the  heat  release  occuring  in  a  region  less  than  about  500  /um  from 
the  tip  ot  the  cavity  in  the  vertical  direction.  The  lateral  extent  of  the  region  is 
about  1500  /xm.  The  maximum  temperature  occurs  at  about  60  fim  from  the 
surf  ace. 

It  was  found  that  values  assumed  for  the  specific  heat  and  the  thermal 
conductivity  of  the  gases  affect  the  solution  significantly  (Ret.  20).  The  effect 
ot  variation  of  the  normalizing  velocity  on  the  solution  is  shown  in  Fig.  7.  The 
figure  dramatizes  the  effect  in  that  the  normalizing  velocity  is  increased  by  a 
factor  ot  3  (cf  Fig.  6,  baseline  case).  The  extent  ot  the  hot  region  increases  and 
the  maximum  temperature  drops.  The  role  of  convective  heat  transter  is  very 
significant. 


1  35 


135 


23 


The  effect  of  the  magnitude  and  the  distribution  of  the  heat  release  was 
also  studied.  For  instance,  it  was  found  that  for  volumetric  heat  release  rates 
above  about  6  x  10^  cal/cm^  sec,  unrealistically  high  temperatures  were 
obtained. 

EXPERIMENTAL  MEASUREMENTS 

The  surface  profile  and  the  interface  burning  rate  have  been  obtained 
experimentally  for  use  in  the  model.  The  burning  sample  is  quenched  and  the 
surface  profile  determined  by  using  an  optical  method.  The  burning  rate  was 
measured  by  using  a  cinephotographic  technique.  The  details  of  these  methods 
are  described  elsewhere  (Ref.  20).  The  measurements  were  performed  on  the 
AP-PMMA  slab  system  burned  at  pressures  of  7.8  and  M.6  atm.  Figure  8  shows 
the  results  of  profile  measurements  for  a  test  at  7.8  atm.  The  vertical  dashed 
lines  in  the  figure  represent  the  conjecture  regarding  the  location  where  the 
walls  might  become  vertical.  The  data  points  which  do  not  follow  this  trend  are 
due  to  relaxation  subsequent  to  softening  and/or  breaking  at  quench.  The 
interface  burning  rate  results  are  shown  in  Fig.  9.  From  the  measurements  of 
burning  rates  during  the  burning  history  of  a  sample,  attainment  of  steady  state 
was  verified. 


MODEL  APPLICATION 

In  this  section,  the  application  of  the  entire  model,  as  described  in  the 
section  "Approach  to  the  Solution,"  is  considered.  The  consideration  is  based  on 
the  use  of  experimentally  determined  surface  profiles  and  burn  rates. 

Let  us  consider  the  combustion  of  an  AP-PMMA  system  at  a  pressure  of  7.8 
atm.  Results  were  computed  for  the  condensed  phase  energy  equation,  gas  phase 
momentum  equation  and  the  gas  phase  energy  equation  based  on  input  parameter 


Figure  8.  Experimentally  determined  surface  profiles 
of  fuel  and  oxidizer  (7.8  atm.). 


values  in  Tables  1,  2  and  3.  The  corresponding  results  are  shown  in  Fig.  10,  11 
and  12.  The  condensed  phase  results  (Fig.  10)  point  to  the  heat  being  transferred 
from  the  oxidizer  to  the  fuel  at  the  vertical  interface.  Figure  11  shows  the 
streamlines  for  the  fluid  solution.  Figure  12  is  a  plot  of  isotherms  for  the  gas 
phase.  Figure  13  shows  the  heat  release  distribution  which  gave  the  above  result. 
The  maximum  temperature  in  the  gas  is  2503°K  and  it  occurs  at  about  1300  nm 
from  the  surface.  All  of  the  heat  release  has  been  assigned  at  locations  below  5 
mm  from  the  surface.  It  may  be  pointed  out  that  these  results  reflect  the 
instance  in  which  overall  convergence  has  been  obtained  between  the  gas 
momentum  and  energy  equations.  The  temperature  gradients  at  the  solid-gas 
interface  have  also  been  matched  satisfactorily.  (See  Ret.  20  for  details.)  The 
total  heat  release  rate  has  been  computed  f  rom  a  thermodynamic  point  of  view 
(Ref.  20).  It  has  been  found  that  about  20%  at  the  heat  release  occurs  at  the 
surface.  The  method  does  not  provide  a  unique  solution,  since  it  is  based  on  an 
iterative  procedure  that  requires  only  approximate  matching  of  modular  solutions 
to  the  conditions  imposed  by  the  input  surface  profile  and  burning  rate. 
However,  several  other  constraints  are  imposed  (Ref.  20)  to  assure  realistic 
results,  and  computer  experiments  have  shown  that  the  solution  (including  heat 
source  distribution)  is  narrowly  constrained  by  the  procedure.  Figure  14  shows 
the  results  for  the  gas  phase  energy  equation  from  a  similar  series  at 
computations  for  a  pressure  of  14.6  atm.  The  maximum  temperature  is  271 1°K 
and  it  occurs  at  a  distance  of  700  /urn  f  rom  the  surf  ace. 

CONCLUSIONS 

A  mathematical  model  has  been  developed  lor  the  combustion  at  the 
interface  of  an  oxidizer-tuel  slab  system.  The  model  is  based  on  two-dimensional 
conservation  equations  and  has  been  solved  numerically.  The  model  uses 


Table  1.  Values  of  input  parameters  for  the  condensed  phase  solution 

(pressure  7.8  atm.) 


PHYSICAL  CONSTANTS 

OXIDIZER 
P|  =  1.95  gm/cm3 

Cj  =  0.3465  cal/gm  °K 

\|  =  0.0009  cal/cm  sec  °K 

5  7 

Ej  =  22000  cal/mole,  Aj  =  3.  x  10  gm/cm  sec 
FUEL 

P2  =  1.185  gm/cm3 

c2  =  0.3  cal/gm  °K 

2  =  0.0005  cal/cm  sec  °K 

E2  =  31000  cal/mole,  A2  =  2.462  x  107 


Slope  of  the  oxidizer  at  the  tip  =  o.O 
Slope  of  the  fuel  at  the  tip  =  0.852 

Interlace  burn  rate  =  0.133  cm/sec  (experimental 
Pressure  =  7.8  atm. 


Table  3.  Values  oi  input  parameters  for  the  gas  phase  energy  solution 

(pressure  7.8  atm.) 


PHYSICAL  CONSTANTS 


\  =  0.0003  cal/cm  sec  °K 

c  =  0.4272  cal/gm  °K 

rO 

Pg  =  0.00092  gm/cm^ 

VQ  =  80.229  cm/sec 
TQ  =  298°K 
hsl  =  -120  cal/gm 
T.  =  2782°K 


(pressure  7.8  atm.;  temp  2782°K) 
p  =  7.8  atm 
r  =  0.133  cm/sec 
h$2  =  314  cal/gm 
Mg  =  26.919 


All  condensed  phase  properties  are  as  in  Table  2. 


Figure  10.  Condensed  phase  isotherm  plot  (pressure  7.8  atm.). 


Figure  11.  Streamline  plot  (pressure  7.8  atm.). 
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experimentally  determined  burning  surface  protiles  and  burning  rate  as  input, 
thereby  avoiding  the  need  tor  modeling  the  actual  t lame  process. 

The  conclusions  from  the  condensed  phase  energy  solution  are  as  follows: 

1.  The  length  scales  over  which  significant  temperature  variations  occur 
have  been  determined. 

2.  The  heat  transfer  occurs  from  the  oxidizer  to  the  fuel  at  the  vertical 
interface  for  the  conditions  studied. 

3.  The  magnitudes  of  the  thermal  conductivities  of  the  solids 
significantly  affect  the  above  conclusions. 

k.  In  general,  the  shape  of  the  surface  profile  affects  the  temperature 
field  more  significantly  than  the  surface  temperature  distribution  (i.e.,  the 
effect  of  surface  decomposition  kinetics  is  secondary).  The  results  point  to 
the  importance  of  obtaining  accurate  experimental  observations  at  the 
profile  at  the  tip  of  the  cavity. 

Solutions  have  been  obtained  for  the  momentum  equations  with  variable 
density  inputs.  The  following  observations  are  made  in  this  context: 

l.  The  inclusion  of  variable  density  leads  to  significant  changes  in  the 
flow  field  (i.e.,  a  constant-density  assumption  leads  to  significant  error  in 
calculation  of  the  flow  field). 

2.  The  coupling  between  the  momentum  and  energy  equations  is  strong. 

3.  Satisfactory  convergence  has  been  obtained  between  these  coupled 
solutions. 

The  gas  phase  energy  equation  has  been  solved.  The  solutions  point  to  the 


following  conclusions: 
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1.  The  solution  is  significantly  affected  by  the  specific  heat  and  the 
thermal  conductivity  of  the  g as  mixture. 

2.  The  role  of  convective  heat  transfer  is  very  significant  in  determining 
the  gas  phase  flame  structure  (i.e.,  representation  of  heat  flow  in  the  gas 
by  a  conduction  term  alone  leads  to  significant  error  in  calculation  of 
temperature  fields  and  heat  source  distribution. 

3.  Realistic  heat  release  rate  distributions  have  been  obtained. 


4.  It  is  believed  that  the  magnitude  and  distribution  of  the  heat  release 
rate  have  been  established  within  certain  limits  by  the  constraints  imposed. 

5.  The  surface  heat  release  is  of  the  order  of  20%  of  the  total  heat 
released  for  the  combustion  of  the  system  at  7.8  atm.  (i.e.,  reasonable  tit 
of  condensed  phase  and  gas  phase  solutions  to  observed  surface  profiles  and 
burning  rates  could  not  be  obtained  without  such  heat  release  at  the 
surface  (even  though  the  tests  were  at  pressures  below  the  self- 
deflagration  limit  of  the  AP.). 
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NOMENCLATURE 


pre-exponential  factor  in  Arrhenius  rate  law  (gm/cm^  sec)  * 


specific  heat  of  solid  phase  (cai/gm°K) 


specific  heat  of  gas  mixture  (cal/gm°K) 


diffusion  coefficient  (cm  /sec) 


activation  energy  (cal/mole) 
volumetric  heat  release  rate  (cai/cm^  sec) 


heat  of  gasification  (cal/gm) 


heat  of  formation  (cal/mole) 


Lewis  number,  X.  /PD  C 


molecular  weight 


unit  normal  vector  at  surface 


nx»  x  and  y  components  of  n 


pressure  (atm.,  other) 


interface  burn  rate  (cm/sec) 


Reynolds  number,  x  v  p/p 


units  in  parentheses  apply  to  corresponding  dimensional  quantities. 


temperature  (°K) 

velocity  in  the  x-direction  (cm/sec) 
velocity  in  the  y-direction  (cm/sec) 
normalizing  velocity  (cm/sec) 
surface  normal  velocity  (cm/sec) 
abcissa 

ordinate 


2 

thermal  dill  usivity  (cm  /sec) 
transformed  y  co-ordinate 
thermal  conductivity  (cal/cm  sec°K) 
viscosity  of  gas  mixture  (gm/cm  sec) 
density  (gm/cm  ) 
shear  stress,  normal  stress 
stream  function 

channel 


gas  mixture 


normal  to  the  surface 


surface 

transition 

wail  (surface) 

x-direction 

y-direction 

conditions  at  infinity 

oxidizer 

fuel 

:s 

dimensional  quantity 

d/dx 

d2/dx2 


vector 
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